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Abstract 

This article is concerned with the asymptotic behaviour, at infinity 
and at the origin, of Green functions of operators of the form Lu — 
— div (AVli) , where A is a periodic, coercive and bounded matrix. 

1 Introduction 

The study of Green's functions for elliptic operators is an important research 
subject. It is linked with many different fields, as for instance homogenization 
[1] [2j [3j [5j [16] , or the study of singular points |T0l [18] . The aim of the present 
article is to provide explicit bounds at infinity for the Green function G of a 
divergence-type elliptic operator with periodic coefficients. Many arguments in 
this paper are already present in the literature in a scattered manner, and our 
main contribution is to put them together in a clear way. Our arguments also 
provide us with explicit bounds on G in the neighbourhood of the origin, where 
G is singular. These latter results are already described in a comprehensive way 
in the literature. 

In all the article, we assume that d > 2 is the dimension of the ambient 
space, and that (here, R dxd is the space of square matrices of size d) the field 



A:R d — > R dxd satisfies 



A is Z d periodic, (1.1) 

A is S — Holder continuous for some 6 > 0, (1-2) 

3a >0, VeeK d , VxeM. d , £ T A(x)£ > a\£\ 2 , (1.3) 

where | ■ | is the Euclidean norm of R d , and 

A g L°° (R d ,R dxd ) . (1.4) 



We want to study the behaviour at infinity of the Green function G associ- 
ated with the operator 

L = — div (AV-) , 
that is, the function G : R d x R d — > K such that 

-div K (A( a; )V x G(a ;! y)) = <5 y ( a; ). (1.5) 

See (|2.ip below for a more precise formulation. By behaviour at infinity, we 
mean the asymptotic of G(x, y) as \x— y\ goes to infinity. This question has been 
widely studied in the literature. According to [TJ Theorem 13] (see also [T5]). 
we have, if d > 3, 

3G, V(i, y ) eM d xM d , \G(x,y)\ <C\x-y\ 2 ~ d . (1.6) 

In addition (see [U Theorem 13]), we have, in the case d = 2, 

3C, V(x,y) et 2 xM 2 , |G(x,y)| < C* (1 + log |x - y|) . (1.7) 

Note that these estimates characterize both the asymptotic behaviour of G at 
infinity (when \x — y\ — > oo) and at the origin (when \x — y\ — > 0). An important 
point here is that many papers consider only the case of Green functions for op- 
erators L defined in a bounded domain (the equation (|1.5[) is then complemented 
by appropriate boundary conditions) . This is the case for instance of [B] and [TU1 
Theorems 1.1 and 3.3]. This is also the case of [TJ Theorem 13], although a re- 
mark following this Theorem indicates that the constant in the estimate can be 
chosen independent of the domain. In [101 Theorem 3.3], bounds are provided 
on G, its gradient and the second derivatives V x V a G, in the case d > 3. A 
remark following that result points out that the constant in the estimate of G 
is independent of the domain, whereas the constants in the estimates of the 
derivatives of G a priori depend on the domain. 

In this article, we also address the question of the decay of the derivatives of 
G at infinity. We have, as proved in Propositions [5] and [7] below (the material 
is present in [TJ, and also in [3]), for any d>2, 

3G>0, V(x,y) eR d xR d , \V x G(x, y)\ + \V y G(x, y)\ < C\x - y^ (1.8) 

and 

3G>0, V(x,y) eR d xM d , \V x V v G(x,y)\<C\x-y\- d . (1.9) 
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A preliminary question, before showing (|1.6p . (|1.7[) , (11.81) and (|1.9p . is the 
existence and uniqueness of G defined by (|1.5p . This question is addressed 
in [Tpl Theorem 1.1], for the Green function in a bounded domain C R d 
with homogeneous Dirichlet boundary conditions. An existence proof is then 
provided for G such that V x G(-,y) G L p (fl \ B r (y)), for any p > d/(d — 1) and 
r > 0. Actually, in [10 , only the case d > 3 is studied, but the existence proof 
carries through to the case d = 2. The uniqueness of G, under the assumption 
that G > 0, is also proved in [TUJ Theorem 1.1] for d > 3. The case d — 2 is not 
covered by their proof. A proof of uniqueness when d = 2 can be found in the 
appendix of [13], both for a bounded domain and for the whole space. 

We finally mention that the case of non-divergence form operators (of parabolic 
and elliptic type) has also been considered, see e.g. [7]. 

The article is organized as follows. In Section [21 we give existence and 
uniqueness theorems for Green functions. In Section [31 we state asymptotic 
properties on G and its derivatives. Finally, we give in Section [4] some remarks 
about possible extensions of the results stated in the present article. 



2 Definition of Green function 

In order to state the existence and uniqueness result for G solution of (11.51) , we 
first write a weak formulation: we look for G : such that 

VyeM d , V^6D(E d ), / (V^(x)f A(x)V x G(x,y) dx = (p(y). (2.1) 



In the sequel, we will need the definition of weak L p spaces, which are special 
cases of Lorentz spaces: for any open subset fl C M. d , for any p € [1, oo], 

LP'™(n) = {/ : Q R, / measurable, < oo} , 

where 



(n) =su P ft M ({xefi, \f( x )\ > t}) 1/p \ 



\L*> 

t>Q 

where \x is the Lebesgue measure. We recall (see e.g. [H p. 8] or [3]) that, for 
any < /3 < p - 1, 

C(p,P,Sl) ||/|U P -^n) < H/IU^O) < ||/|Uo(fi), (2.2) 



with c(p, p, o) = m fc(n)) hf^t . 

Theorem 1 (Existence and uniqueness of G, d > 3). Lei > 3, and assume 
that A satisfies (|1.3p and p.4p . Then, equation (|2.ip /ias a unique solution in 

L™ (R d ,W*ll c {®. d j) such that 

lim G(x,j/)=0. (2.3) 

|a: — v|— >-oo 
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Moreover, G satisfies the following estimate: 

Vg<^4i, VyeM^, G(;y)eW^(R d )nW^(R d \{y}) (2.4) 

and 

3C, V{x,y) eR d xR d , < G(ir, y) < C\x - y\ 2 ~ d . (2.5) 

Proof. First, note that, according to [HI Theorem 8.24], the function G is Holder 
continuous with respect to x and y whenever x ^= y. The same property holds 
for Gr defined below. 

Let R > 0. We first define Gr as the Green function of the operator 
— div(AV-) on the ball Br — Br(0) with homogeneous Dirichlet boundary 
conditions, that is, 

VyeB R , VipeV(BR), f (V<p{x)f A(x)V x G R {x,y)dx = <p(y), (2.6) 

Jb r 

and Gii(x,y) — if \x\ — R. Applying (TUJ Theorem 1.1], we know that such a 
Gr exists, and satisfies 

VyeB R , \\G R (.,v)\\ * <C, (2.7) 

L*- 2 ' (Br) 

V?y G B R , \\V x GR(;y)\\ < C, (2.8) 

i 1 ^ 1 (Br) 

and 

V(x, y)eB R x B R , < G R (x, y) < ^™ (2.9) 

\x — y\ 

where C > does not depend on R and y. 

Next, we note that if R' > R, then, due to the maximum principle, we have 
Gri > Gr in Br x Br. Thus, Gr is a non-decreasing function of R. With the 
help of ()2.9|) . this implies that the function Gr converges almost everywhere to 
some function G, defined on ~R d x ~R d , and that satisfies (|2.5[) . This implies (|2.3[) . 
In addition, we deduce from flJjU) that G fl converges to G in £f oc (K d x R d ), for 
any p < d/(d — 2), and that, for any y £ K d , the function Gr(-, y) converges to 
G{-,y) in Lf oc (R d ), for any p < d/(d - 2). 

In view of (|2.8j) and (|2.2|) . we see that, for any bounded domain il C M. d , 
and for any q < d/(d — 1), there exists G(f2, g, d) such that 

Vi? s.t. n C B R , Vy G B Bl ||V B Gji(.,y)|| £ , (n) < G(fi, <?,<*). 

Hence, extracting a subsequence if necessary, V x Gfl (•,?/) converges weakly in 
(L 9 (f2)) d to some T £ (L 9 (r2)) d . Recall now that Gr(-, y) converges to G(-, y) in 
Lf oc (M d ), for any p < d/{d — 2). Hence T = V x G\n, and V x G_R(-,y) converges 
to V^G weakly in (L q (Q)) d , for any bounded domain Q and any q < d/(<i — 1). 
Passing to the limit in ()2.6[) . we see that G is a solution to (|2.1[) . 

Finally, the bounds (|2~7|) and (|2TS|) imply, together with (j2T2"j) . that G £ 

W^,k>c( Rd )) • We have thus Proved the existence of G. 
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Property (|2.4[) is proved in [TU1 Theorem 1.1], and its proof does not depend 
on the fact that the domain used there is bounded. Note that we have already 
proved part of this property. Indeed, as pointed above, for any y G R'', we have 
G(-,y) 6 Ll c (R d ) for any p < d/(d - 2) and V x G(-,y) € (L q loc (R d )) d for any 
q < d/(d - 1), thus G(; y) G W^{R d ) for any q < d/(d - 1). 

In order to prove uniqueness, we assume that G\ and Gi are two solutions, 
and point out that H = G\ — G2 satisfies div^ (AV X H) = for any y € R d . 
Fixing y, we apply the corollary of [20l Theorem 4], which implies that, if H 
is not constant, then sup {H(x, y), \x — y\ = r} — inf {H(x, y), \x — y\ = r} must 
grow at least like a positive power of r as r —> 00. This latter behaviour is in 
contradiction with (|2.3I) . Thus H = G\ — G2 is constant, and (|2.3p implies that 
G\ = G'i- 

Note finally that the corollary of [20l Theorem 4] is stated in the case when 
A is symmetric, but the same result holds in the non-symmetric case. Indeed, 
Harnack's inequality is still valid in such a case, see e.g. [SJ Theorem 8.20], 
Theorem 5.3.2] or [T3]. □ 

Theorem 2 (Existence and uniqueness of G, d — 2). Let d — 2, and assume 
that A satisfies (|1.3[) and (|1.4p . Then, equation (|2.1[) /ias a unique (up to the 

addition of a constant) solution in L^ loc MR rf , W^'j^fR^H smc/i i/iai 

3C*>0, V(s,|/) S M d x K d , |G(x,y)| < C(l+ |log|ar-»||) . (2.10) 

Moreover, G satisfies the following estimate: 

V-z<2, \fyeR d , G(;y)eW^(R d )nW^(R d \{y}). (2.11) 

Proof. The proof of this result may be found in the appendix of |14j . However, 
for the sake of completeness, we provide an alternative proof. This proof, in 
contrast to that of [M], relies on basic tools of analysis of PDEs. 

We use the same strategy as in the proof of Theorem [1] defining first the 
Green function Gr of the operator L on Br. However, we cannot simply apply 
the results of [10] to define Gr, as those results are stated in dimension d > 3. 
It is possible to adapt the proof of [TUJ Theorem 1.1] to the two-dimensional 
case, but a simpler proof consists in following the approach of Section 6 of [6] . 
These results give the existence and uniqueness of Gr solution to (I2.6P in the 
ball B R = B R (0) C JR 2 , with the homo geneous Dirichlet boundary conditions 
GR(x,y) = if \x\ = R, in W 1,p (Br) for any p < 2. In addition, it is shown 
in [51 Section 6] that estimate (12.81) holds, namely 

VyeB R , \\V x GR(;y)\\ L 2,~ (Ba} <C (2.12) 

for a constant C independent of R and y. 
Step 1: passing to the limit R — >• 00 on Gr, 
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Consider the domain = Br>, with R' fixed, and consider next R > R'. 
Applying (|2.2[) to \7 x Gr(-,h) on Q, we see that (|2. 12|) implies that V x Gr is 
bounded in (L q (BR> x -Br/)) 2 for any q < 2, independently of i?. Hence, ex- 
tracting a subsequence if necessary, V x Gr converges weakly in (L q (Bjj/ xBr')) 2 
to T 6 (L q (Bni x Bui)) 2 . Now, we have, in the sense of distribution, 

d Xl d X2 GR = d X2 d Xl Gn. 

This property passes to the limit, so that d Xl T2 = d X2 T\. This implies that 
T = V X G for some G e VF 1 ' 9 {B R , x £r/). Next, we point out that this limit 
does not depend on R 1 in the sense that if R" > R' , then V X G' obtained 
in Bri is equal to V X G',' B ; , where V X G" is obtained in Hence G £ 

Wloc ( R2 xR2 ) c L 9,ioc (^f^xicO* 2 ))- Passin g to the lim it in we 
obtain that G is a solution to (12.11) . Until now, the function G(-,y) is only 
determined up to a constant. We fix this constant by choosing G(-, y) such that 

/ G{x,y)dx = 0. (2.13) 
Jb 1 { v ) 



To prove the existence of a function G satisfying the claimed properties, we 
are now left with showing that the function G that we have built satisfies (|2.10p 
and (|2~TT|) . 

Step2: proving that G satisfies (12. lip 

By construction, we have G(-,y) € W 1,9 (f2), for any q < 2 and any bounded 
domain fl The proof of the fact that G(-,y) £ W 7 ^ 2 (M 2 \ y) follows the same 
lines as the proof given in [TUl Theorem 1.1], which does not depend on the fact 
that the domain used there is bounded, nor on the fact that the dimension there 
is d > 3. We thus have proved (|2 . 1 1 [) . 

Step 3: proving that G satisfies (|2.10p 

We first infer from (|2.12l) and (I2.2p that, for any bounded domain ft C Br 
and any y £ Br, we have 

\\V x G R (;y)\\ LHn) <C 



for a constant G independent of R, f2 and y. Since V x Gr(-, y) weakly converges 
to V x G(-,y), we deduce that 

|V B G(.,y)|| i i (n j<C (2.14) 



for a constant G independent of fl and y. Note that this implies that G G 
L ^ioc ( Rd > W^ioc^)) > as claimed in the theorem. 



G 



Second, we apply Poincare-Wirtinger inequality to G(-,y) on the set Bi{y): 
using (|2.13p . we have 



G(x,y)\dx<C \V x G{x,y)\dx. 

Bl(y) JB 1 (y) 

Applying (|2.14[) with VL — B\(y), we deduce that 

\G(x,y)\dx < C, (2.15) 

Bi(y) 

where C does not depend on y. 

We next define, for any R > 0, the function 



/(^) = ^rs/ \G(x,y)\dx, 

>dB R (y) 



2ttR 



where dx denotes the Lebesgue measure on the circle dBf>(y). Note that / 
depends on y, but we keep this dependency implicit in our notation. In the 
sequel of the proof, we first show a bound on / (step 3a) , and then deduce from 
that bound a bound on G (step 3b). 

Step 3a: bound on / 

We have, for any R > R' > 0: 



\f(R)-f(R')\< j R \f{r)\dr< [ R ^- j |V a 
Jr> Jr> 2nr J dBr ( v ) 



G{x, y) \ dxdr 



1 f \/R 2 - R' 2 R 
< ^7 / \V x G(x,y)\dx < C^— = C-, (2.16) 

JB R (y)\B Rl {y) « « 

where we have again used (|2.14p and where the constant C does not depend on 
y. This implies that f(R) is bounded independently of R and y for R G (1/2, 1). 
Indeed, for such an i?, we rewrite (gHH) as /(i?) < f{R') + CR/R' (recall that / 
is non- negative) , and integrate with respect to R' between 1/4 and 1/2, finding 



-J{R) < [ f(R')dR' + CR. 

4 Jl/A 

Using (|2.15p . we infer 



f(R) < C, (2.17) 



for some constant C independent of R and y. Next, we consider two different 
cases: R > 1 and R < 1/2. 

• Case R > 1: in such a case, we define p€ N such that 

1 i? 

- < — < 1, 

2 2p ~ 
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that is, p is the integer part of which reads < p < + 1. 

We then apply (pHo]) with R = 2~iR, R' = 2~i- 1 R, finding 



2J 7 ' \ 2i 



where C is a constant which does not depend on R, j, nor on y. We sum 
up all these inequalities for < j < p — 1, and obtain 

/(fl)</(^)+Cp. 

Recalling (|2 . 1 7[) and the definition of p, we infer 

\f(R)\ <C (1 + | log(i?)|), (2.18) 
where C is independent of R and y. 

• Case R < 1/2: the approach is similar to the preceding case. We define 
p EN such that 

\ < 2 p i? < 1, 

that is, p is the integer part of -^f - 1. We apply ([2~T6]) with R' = VR 
and R = 2^ +1 R, finding 

\}{VR) - f{V +1 R)\ < C. 

We sum this with respect to < j < p — 1, and find that (|2.18j) is again 
valid in this case. 

Collecting the result of the above two cases, we find that 

Vi?>0, |/(iZ)|<tf(l + |log(iJ)|), (2.19) 

where the constant C does not depend on R nor on y. 

Step 3b: bound on G 

We first make use of (|2.19| to obtain a bound on the L 1 norm of G in any 
annulus. For any /3 < 7, we indeed have 

rl 

\\G(;y)hi(B^(y)\Be(y)) =2ir rf(r)dr, 
hence, using (|2.19|) . we obtain 

\\G(;y)\\ Ll{B ^ y) \ Blj{y)) <C T r(l + \\og(r)\)dr. (2.20) 
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Consider now R>l/2. Then 3R>2R> 1, and (|2~2T)1) implies 

r3R 

Vi?>l/2, \\G(;y)\\ L i {B3RmB2R[v)) <C r(l + log(r))dr 

J2R 

< 3CR 2 (1 + log(3i?)) < CR 2 (1 + \log(R)\) , (2.21) 
for some C independent of R and y. In turn, if R < 1/3, then (I2.20P implies 

ViZ<l/3, \\G(;y)\\ L , {B3R{y) \ BMy)) <C / r(l-log(r))dr 

< 3C*i? 2 (1 - log(3i?)) < CR 2 (1 + |log(i?)|) , (2.22) 

for some C independent of R and y. 

Next, we recall that, according to Sobolev imbeddings (see for instance [HI 
Theorem 7.10]), we have 

Vp<2, VueW hp (R 2 ), H| je. < C p \\Vu\\ LHm . 

L 2 -p(R 2 ) v ' 

We apply this inequality tou = G(-, y)xRi where \R 1S a cut-off function satis- 
fying 

X r eV(R 2 ), |Vxh|<-, XR = outside B 3R (y), X R = l™B 2R (y). 
We find, for p = 1, that 

\\ G {->y)\\L2(B 2R {y)\B R (y)) < IMU 2 (K 2 ) < C|| Vu[|£l( H a) 

c 

< C\\V x G(-,y)\\ L i(B 3R (y)) + -^\\ G (^y)\\L^B 3R ( y )\B 2R (y))- ( 2 - 23 ) 

The first term of the right-hand side is bounded using (|2.14l) , that yields 

\\\7 x G(- 7 y)\\ LHB3R{y)) <CR. (2.24) 

The second term is bounded using (ET2"0"|) . ([2T2"Tj) and ([2T2"2"j) . If i? > 1/2 or 
i? < 1/3, we indeed see from (pT2Tj) and (|2~2"2l that 

^l|G(-,l/)|Ui( BM („)\ Bail ( v )) < Cii (l+|log(iJ)|). (2.25) 

In turn, if 1/3 < R < 1/2, then we deduce from ((11201) that 
C 

-ft\\G{-,y)\\mB 3R (y)\B 2R (y)) < 3C 1 1 G(- , ?/) || L i { b 3/2 (y)\B 2/3 (y)) < C, 

and hence (|2.25p is again valid. 

Collecting (|2T23|) . (|2~24l) and (|2~25]) . we obtain 

\\ G \\Li(B 2R (y)\B R (y)) < <?-R + CR |k>g(i?)| . 



Finally, we apply Theorem 2] (see also [HI Theorem 8.15]), which implies 
that, for any v € W lo ' c (M 2 ) such that Lv = in B iR {y) \ B R / 2 {y), we have 

C 

sup v < -\\v\\ L 2 {B2r{v) \ Br{v)) . 

B 2R (y)\B R (y) H 

Applying this to G(-,y) and —G(-,y), we find 

sup |G(.,i/)|<C(l + |log(iJ)|). (2.26) 

B 2 R(y)\B„(y) 

The function G hence satisfies (|2.10l) . This concludes the proof of the existence 
of a function G satisfying the properties claimed in Theorem [2j 

To prove the uniqueness of G (up to a constant), we follow the same argument 
as in the case d > 3 (see Theorem[JJ). Assume that G\ and G2 are two solutions. 
We point out that H = G\ — G2 satisfies div^ (AV X H) = for any y G K 2 . 
Fixing y, we apply the corollary of |20[ Theorem 4], which implies that, if H 
is not constant, then sup {H(x, y), \x — y\ = r} — inf {H(x, y), \x — y\ = r} must 
grow at least like a positive power of r as r — > 00. This latter behaviour is in 
contradiction with (|2.10p . Thus H = G\ — G2 is constant. This concludes the 
proof of Theorem [2] □ 

Remark 3. The above proof can be adapted to the case of the Green function 
Gu on the bounded domain Br, i.e. the solution to (|2.6p . We hence obtain 

V(x,y) £B R x B R , \G R {x,y)\ < C R + C|log \x - y\\, 

thus recovering the result o/|3 Section 6]. Note that the constant Cr in the above 
bound a priori depends on R. Think indeed for instance of the case L = —A, 
where Gr(x, 0) = — log + log R. 

3 Asymptotic behaviour 

We now give some results about the asymptotic behaviour (at infinity and at the 
origin) of the Green function G. First, we note that, collecting (12.51) and (|2.10[) , 
we have the following: 

Proposition 4. Assume that A satisfies (|1.3p and (|1.4|) . Then, the Green 
function G of the operator — div(AV-) (namely the solution to (|2.ip ) satisfies 

3C- > 0, V(,, S ) e R* x W, \G(* M < " "» J * (3.1) 

As we pointed out in the introduction, this result is well-known for bounded 
domains [TJ [6j [101 [15l [18] . However, almost all results are limited to this case, 
except for [TJ Theorem 13], for which "in spirit", the domain is infinite due 
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to the scaling with respect to e — > 0. The articles [T5] and [T5J Section 10] 
also consider the case of unbounded domains (see also a remark following [Tol 
Theorem 3.3]), but do not consider the case d = 2. Finally, the appendix of Q3] 
treats the case of M 2 . 

Next, we give results on the gradient of G. 

Proposition 5. Assume that A satisfies (|1.1[) . (|1.2[) . (|1.3p and p.4[) . Then 
the Green function G associated with L = — div (AV-) satisfies the following 
estimates: 

3C>0, VxeR d , VyeR d , \W x G(x,y)\ < (3.2) 

\x-y\ a 

30 0. VxeR d , Vy &R d , \V y G(x,y)\ < -. -. (3.3) 

\x — y\ a 

Similar results are given in |10| Theorem 3.3], in the case of bounded do- 
mains. 

Proof. We start with the case d > 3, and apply [U Lemma 16] to G as a function 
of x, which implies that 

VxeM d , VyeIR'', Vr<\x-y\, 

\\V*G(',v)\\L-( Br/a (x)) < - \\G(-,y)\\ L ~(B r ( x )), (3.4) 
where C depends only on ||A|| c o,«, 5, a and d. Using (|3.ip . we thus obtain 

\V x G(x,y)\<- sup 1 — r ^. (3.5) 

r zeB r (x) \ z - y\ 

Note that we have used \V x G(x,y)\ < \\y x G(-,y)\\[ J ^(B r/2 (x)), which is true 
only almost everywhere. However, changing the function on a set of measure 
zero if necessary, it is possible to assume that this inequality holds everywhere. 
Taking r = \x — y\/2, we have, for any z e B r (x), 

\x-y\ < \x-z\ + \z-y\ < r + \z - y\ = -\x -y\ + \z-y\. 



We hence deduce from (13.51) that 



\VxG(x,y)\ < - 



C ( 2 \ d ~ 2 2 d - 1 C 



r \ \x — y\j \x — y\ d 1 



This proves (pT2|) . 

Next, in order to prove (|3.3|1 . we point out that G*(x,y) :~ G(y,x) is the 
Green function of the operator L* defined by 

L*u = - div (A T Vu) . (3.6) 
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A proof of this fact0 is given in [TU1 Theorem 1.3] and [BJ Theorem 1] in the 
case d > 3, and this proof carries over to the case d — 2. Hence, applying (I3.2[) 
to G* we deduce 



We turn to the case d = 2. The estimate (13.41) is not sufficient here, since 
G(x, y) is not bounded as \x — y\ — > oo. Instead, we use the same trick as in the 
proof of [TJ Theorem 13], using (|3.2[) for d — 3. For this purpose, we introduce 
the operator L defined on lf 1 (R 3 ) by 



Lu = — div x (A(x)V x u) — <9 t 2 u, 



(3.7) 



where iGl 2 and t S R. Let G be the associated Green function. According to 
the above proof and to (|3. 1[) , we have 



G(x,t,y,s) 



< 



C 



and 



V x G(x,t,y,s) 
Next, we set, for any x and y in 

G K {x,y) = 

We deduce from (j3~gf that 

\V x G K (x,y)\ <C 



x-y\ + \t- s 
d t G(x,t,y,s) < 



C 



x - y\ 2 + (t- s) s 
, with x y, 



(3.8) 



G{x,t,y,0)dt. 



dt 



Cir 



y\ 



(3.9) 



for a constant C independent of k, x and y. Hence, V X G K is bounded in 



uniformly with respect to n, for any p < 2. Thus, for any R > 0, 



extracting a subsequence if necessary, V X G K converges weakly in (L p (B rX B p)) 2 
to some T S (L p (Br x Br)) 2 . Now, we have, in the sense of distribution, 



^x\^x^G K d X2 d xl G K . 



This property passes to the limit, so that d xl T 2 — d X2 T\. This implies that 
T = V X G for some G G ^^^(Bfl x B R ). We next point out that the limit G 
does not depend on R, in the sense that if R' > R, then G defined on Br> x Bri 



as above satisfies \7 X G = S7 X G. We thus have GeL 



oo 

y ,loc 



1,1 

x,loc 



Note also that p.9[) implies that, for any y G M 2 , the function V a; G' K (-,y) 
is bounded in L^ oc (R 2 ), uniformly with respect to k, for any p < 2. Thus, 



lr The main idea of the proof consists in choosing the test function <p(x) = G(z,x) in l|2.1|l . 
for any z £ R d , and next multiplying 1 12.11 1 by an arbitrary function f(y) and integrating over 
y. However, as the function G(z, ■) does not belong to X>(R d ), some regularization arguments 
are in order. 
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for any bounded domain Br, extracting a subsequence if necessary, V X G K (-, y) 
converges weakly in {L p {Br)) 2 , and, by uniqueness, V x G K (-,y) converges to 
V x G(-,y) weakly in {U>{B R )) 2 . 

At this point, G(-, y) is only determined up to an additive constant. We now 
fix this constant (and hence uniquely defined G (•,?/)) by assuming that 

G(x, y) dx = 0. 

Bi{y) 

In the sequel, we show that G satisfies all the properties of Theorem [5] By 
uniqueness of the Green function G up to an additive constant, we will obtain 
that G = G up to a constant. We will then deduce bounds on VG from the 
bounds we have on VG. 

We first show that G satisfies ([27T]) . Consider ip G D(R 2 ) and tp G X»(K). 
Considering the test function ip(t)(p(x) in (|2.1j) . we see that the Green function 
G satisfies the weak formulation 

tp(t) {V<p{x)) T A(x)W x G(x, t, y, 0) dxdt 

<p(x) V>'(i) 9 t G(ar, t, y, 0) dxdt = (p(y)i/>(0). 

Consider tp such that i(j(t) = 1 whenever \t\ < k, ip(t) = whenever \t\ > 1 + k, 
and maxdlV'IUoo(M), IIV>'IU~(R)) < L We have 

(V<p{x)) T A(x)V x G K (x, y) dx + e x («) + ea(/t) = <p(y), (3.10) 

with 

ei(«) = [ [ i>{t){V<p{x)) T A{x)\7 x G(x,t,y,0)dxdt, 

JR 2 JK<\t\<l+K 



e 2 («) = / / <p(x)ip'(t)d t G(x,t,y,0)dxdt. 

JR 2 J K <\t\<l+K 

Let us now bound from above e\ and e 2 - Using (|3.8p . and introducing a compact 
A" C R 2 containing the support of tp, we have 



|ei(«)| < ||V|U°o||A|Uoo||V^|U=o / / V x G{x,t,y,0) 

JK Jk<|*|<1+k 

< ||^|| £ -||A|| £ -||V^||l- f [ , -jj— J dard* 

•As: ./«<|t|<i+« 1^ — 2/1 + * 



dxdt 



t\<l+n 
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Hence, ei(fc) vanishes when k — > oo. Likewise, e2(«) also vanishes when k — > 
oo. Passing to the limit k — > oo in p.lOp . and using that V x G K (-,y) weakly 
converges to V x G(-,y), we deduce that, for any <p e £>(R 2 ), we have 

/ (V^(x)f A(x)V x G(x, y) dx = <p(y). 

We have thus obtained that the function G e L£° loc (r 2 , (M 2 )) satis- 

fies (|2.ip . Assume now that G also satisfies (12 . 10[) . Then, according to the 
uniqueness of G (see Theorem [2]), we have V X G = \7 X G. 
In turn, we deduce from p.9[) that 

\V x G(x,y)\ = \V x G(^v)\<T^ L -r (3-H) 

This hence proves the estimate p.2p in the case <i = 2. 

To prove p.3[) in the case d = 2, we again use the fact that G(y, x) is the 
Green function of L* defined by (|3.6I) , so the estimate (|3.2[) that we have just 
shown implies (13.31) . 

There only remains to prove that G satisfies (j2.10[) . To this end, we note 
that p. lip implies the estimate (I2.14p . for fi a ball or an annulus of the form 
B2R \ Br. Hence, the end of the proof of Theorem [2] applies here, leading 
from (j2~14l) to (pT26|) . which implies that G satisfies ([2~10| . □ 

Remark 6. The above arguments indicate two different proofs for the existence 
of G in dimension two: the first one consists in defining the Green function on 
the bounded domain Br, and then letting R — > 00 ; as it is done in the proof of 
Theorem [H The second strategy uses the three-dimensional Green function G 
of the operator L defined by p.7p . One integrates G with respect to the third 
variable, finding a Green function for the operator L in dimension two. This 
approach is used in the proof of Proposition [31 

Note also that Proposition^ is proved under stronger assumptions than The- 
orem^ 

We next prove upper bounds on V x \7 y G. 

Proposition 7. Assume that A satisfies (TTTTI) . (TO)) . (TO)) and (fi~4"l) . Then 
the Green function G associated with L = — div (^4V-) satisfies the following 
estimate: 

3G>0, VxeW 1 , VyeR d , \V x V y G(x,y)\ < C (3.12) 

\x - y\ d 

Here again, similar results for the Green function in a bounded domain are 
given in the literature, for instance in jTOJ Theorem 3.3]. 

Proof. We have, in the sense of distribution, 

- div x (A(x)V x V y G(x, y)) =0 in B s (y) c , for any S > 0. 
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We can thus apply [TJ Lemma 16], and obtain, as in (|3.4p . that 



VxeR d , VyeM d , W<\x-y\, 

C 

||V a; Vi,G(-,y)||z,oo( Bj . /2 ( a .)) < —\\VyG{-,y)\\ L ^ {Br{x)) . 

Using ([O]) . we deduce (ETT2)) . □ 

Using arguments similar to those used to prove Propositions[5]and[71 we also 
show the following result on the Green function Gr of the operator — div (AV-) 
on the bounded domain Br with homogeneous Dirichlet boundary conditions. 
The interest of this result is the independence of the obtained bounds with 
respect to the size of the domain Br. 

Proposition 8. Assume that A satisfies (jl.ip . ([L2]) . (fl~3|) and (fL~4"l) . Let G R 

be the Green junction of the operator — div (^4V-) on Br with homogeneous 
Dirichlet boundary conditions (namely, Gr is the unique solution to (12 .6[) with 
the boundary condition GR(x,y) — if \x\ = R). 

Then, there exists a constant C such that, for any R > 0, 

V(x,y)eB R xB R , \V x G R (x,y)\ < ^ ZT , (3.13) 

\x-y\ d 

V(x,y)£B R xB R , \V y G R {x,y)\<- (3.14) 

\x - y^- 1 

V(x,y)GB R xB R , \W x V v G R {x,y)\ < j^r^p - (3-15) 



4 Extensions 

First, it should be noted that, assuming further regularity on the coefficients 
of the matrix A, it is possible to prove more precise decaying properties of the 
Green function. This was proved in |21) . 

Next, it is clearly possible to adapt the technique of [B] and [5] (see also [TTJ 
12J) to treat the case of systems of elliptic PDEs. This case is also considered 
in [J. 

Another question is the extension of the present results to the case of non- 
periodic coefficients. This is, for instance, what is done in [TU] and [5], in the 
case of bounded domains. However, some of the estimates we have used here 
(in particular (|3.4|l ) do rely on the fact that the matrix is periodic. Thus, the 
extension is not straightforward. 

Finally, it should be possible to extend our results to the case of piecewise 
Holder coefficients. For instance, gradient estimates for elliptic equations with 
such discontinuous coefficients are derived in |17) . It is probably possible to use 
them in the setting of the current article, but such a work remains to be done. 
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